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Faraday ghosts: depolarization canals in the Galactic radio 
emission 
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CLf 1 INTRODUCTION 



ABSTRACT 

Narrow, elongated regions of very low polarized intensity - so-called canals - have re- 
cently been observed by several authors at decimeter wavelengths in various directions 
in the Milky Way, but their origin remains enigmatic. We show that the canals arise 
from depolarization by differential Faraday rotation in the interstellar medium and 
that they represent level lines of Faraday rotation measure RM, a random function 
of position in the sky. Statistical properties of the separation of canals depend on the 
autocorrelation function of RM, and so provide a useful tool for studies of interstellar 
turbulence. 
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Radio polarization maps of the Milky Way often exhibit a 
network of apparently randomly oriented, narrow (less than 
the beam width), elongated regions where polarized inten- 
sity vanishes or is greatly reduced. These structures, now 
known as 'canals', have been reported by a number of au- 
thors and belong to the class of structures visible in po- 
larized emission that do not have counterparts in the total 
radio intensity (Duncan et al. 1997; Gray et al. 1998, 1999; 
Uyaniker et al. 1999a; Haverkorn, Katgert & de Bruyn 2000; 
Gaensler et al. 2001). In Fig. [I] we show an example of canals 
in the field of the nearby galaxy M31. Some of the canals 
are nearly closed, forming cells, or at least have pronounced 
curvature, and within the cells the polarization angle often 
hardly varies. Other canals are more or less straight, and 
often occur in nearly parallel pairs. The appearance of the 
canals depends on the wavelength of observation (Haverkorn 
et al. 2000). An example of similar features was discussed 
by Scheuer, Hannay & Hargrave 1197711 . 

The canals were first identified as distinct features in 
polarization maps by Duncan et al. 1199711 and Uyaniker et 
al. <!1999all . The latter authors propose that the canals may 
be caused by filaments in thermal plasma and/or suitable 
magnetic structures. As shown by Haverkorn et al. {2000), 
the polarization angle of the radio emission flips by about 
90° across a canal, which prompted these authors to sug- 
gest that the canals arise from beam depolarization due to a 
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discontinuous distribution of the foreground Faraday rota- 
tion measure RM caused by abrupt changes in the magnetic 
field direction. Gaensler at al. 1200111 support the idea of 
beam depolarization, but note that the foreground RM does 
not exhibit abrupt changes across the canals, and therefore 
attribute these structures to a highly nonuniform magnetic 
field in the synchrotron source. It seems, however, highly im- 
plausible that the distribution of magnetic field should be 
discontinuous to that extent, and that either the variations 
of RM have just the right amplitude to rotate the polar- 
ization angle preferentially by 90° or the magnetic field di- 
rection changes preferentially by right angles. Furthermore, 
neither of these interpretations can easily explain why the 
pattern of canals changes with wavelength. 

The canals apparent in the field of M31 (see Fig.0 must 
originate in the Milky Way because their network extends 
far beyond the image of M31. An eyeball estimate of the 
typical separation of the canals is about 5', corresponding to 
a linear scale of about 1 pc at a distance of 1 kpc. Figure ^ 
confirms that the polarization angles flip by about 90° across 
the canals. 

Uyaniker et al. 119981 I1999at mention the possibility 
that the canals may be due to Faraday depolarization with- 
out giving any details. We attribute the canals to depolar- 
ization along the line of sight rather than across the beam. 
This was first proposed by Beck ( 1999] I2001t who showed 
that the canals can arise from depolarization by differen- 
tial Faraday rotation, and so do not require any discontinu- 
ities or even strong gradients in the parameters of the in- 
terstellar medium (ISM). In this letter we substantiate this 
idea and argue that the canals represent level lines of RM 
where RM = RMo = n7t/(2A 2 ) with A the wavelength and 
n — 1, 2, 3, ... . Then we discuss how statistical properties of 
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Figure 1. (a) Polarized intensity at A20.5cm in the field of the galaxy M31, obtained by smoothing the high-resolution map of Beck, 
Berkhuijsen & Hoernes (19983 to an angular resolution of 3' (no corrections for missed spacings were made) . The frame measures 130' X 60' 
with the major axis of M31 oriented horizontally, and is centred on the nucleus at Galactic coordinates (I, b) ?s (120°, —20°). Because of 
strong depolarization, polarized emission from M31 is only visible as slight enhancements of emission along a large, elliptical ring and 
near the nucleus, better visible in total intensity in (b). Canals are visible as an irregular network of narrow, low-intensity regions (light 
shades of grey). As this network extends to the borders of the map, i.e. to regions outside of M31, the canals and much of the polarized 
emission must originate in the Milky Way. (b) The distribution of the B-vectors of polarized emission (observed E-vectors rotated by 
90°) at A20.5cm in the field of M31 shown in (a). The angular resolution of 3' is indicated in the lower left corner. The length of the 
vectors is proportional to the polarized intensity with 3' corresponding to 2 mJy/beam area. The canals seen in (a) are visible as narrow 
lanes (less than 3' wide) that are free of vectors. Note the sudden changes in the polarization angle by about 90° across the canals. 
Contour levels of 5, 10, 15, 20, 30, 40 and 60 mJy/beam area show the total power emission from M31 at A20.5cm. 



the ISM can be determined from the angular separation of 
the canals. 



2 DIFFERENTIAL FARADAY ROTATION 

Differential Faraday rotation is one of the simplest and best 
known depolarization mechanisms which arises where rela- 



tivistic and thermal electrons fill the same region of mag- 
netized medium, so that synchrotron emission and Faraday 
rotation occur together. The degree of polarization p then 
varies with the observable Faraday rotation measure RM 
and A as (Burn 1966, Sokoloff et al. 1998) 
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Figure 2. Differential Faraday rotation under a range of intrin- 
sic Faraday rotation measures. Each horizontal strip represents 
one telescope beam with the observer to the right; the continuous 
curve represents the function 2A 2 RM(a;) that varies across the 
sky (along x); the intrinsic polarization angle is uniform (vertical 
double arrows on the left) . Emission generated on the x-axis can- 
cels the emission generated in a layer where the former is rotated 
by 7r/2, i.e., where 2A 2 RM(x) = 7r/2. Regions whose emissions 
mutually cancel are indicated by perpendicularly hatched boxes. 
Only the emission from unshaded regions reaches the observer. 
The arrows on the right show the observed polarization plane 
which has been Faraday rotated in layers to the right of the vis- 
ible (unshaded) region to give the polarization angle i/> in the 
range indicated. Note that layers to the right of the source act as 
a Faraday screen, and so the rotation angle is 2A 2 RM. The canals 
occur where the polarized emission is fully cancelled (regions la- 
belled B). The emission observed at values of 2A 2 RM(a:) slightly 
smaller than nir (regions A) and slightly larger than rnr (regions 
C) originate in different layers. The difference in 2A 2 RM(x) be- 
tween regions A and C is 7r/2. Therefore, emission from region A 
is rotated by it/2 before reaching the observer, whereas emission 
from region C is rotated only slightly. Thus, the polarization an- 
gles observed on the two sides of a canal differ by 7r/2. A canal 
that occurs when 2A 2 RM(x) = n-K with n = 1 is shown near the 
top, and with n = 2, near the bottom. 

and so vanishes when ]RM| = RMo, where 

2RM A 2 = nn , n = 1, 2, . . . , (2) 

with po ~ 0.75 the intrinsic degree of polarization This 
expression is applicable to a uniform slab and to a few 
more complicated configurations discussed by Sokoloff et al. 
( 1998 ), and is sufficient for our purposes here. It is important 
to note that the angle of polarization suffers jumps by 7r/2 
at the same values of RM where p vanishes [see, e.g., Fig. 2b 
in Burn flM)l and Sect. 3.1 in Sokoloff et al. {1558^ . The 
nature of depolarization by differential Faraday rotation and 
the origin of the jumps are further illustrated in Fig. [5] 

If RM varies across the sky, polarized intensity will 
vanish at those positions x in the sky plane where 
2|RM(cc) | A 2 — nn, i.e., for a given wavelength, at the criti- 
cal level lines of RM where |RM(x)| = RMo. For example, 
RMo fa 37nradm -2 for A = 20.5 cm. The loci of vanish- 
ing polarization are geometric lines, i.e., they are infinitely 
thin, hence always narrower than the beam, as are the canals 



(Gaensler et al. 2001; see Fig.^). The position of the level 
lines changes with wavelength, as does the network of canals. 
Extended polarized emissions seen on different sides of a 
canal come from different regions along the line of sight (see 
Fig-HJ- Hence, the polarization angles differ by 7r/2 as RM 
crosses the critical level RMo, i.e., across the canal. Faraday 
depolarization does not affect the total synchrotron emis- 
sion, so canals are not visible in total intensity. We believe 
that these arguments are convincing enough to attribute the 
canals to depolarization by differential Faraday rotation. As 
they are not directly related to any structures in the ISM, 
we call them 'Faraday ghosts'. 

We now recall that in the turbulent interstellar medium 
RM is a random function of position. As discussed by War- 
die & Kronberg 11974t and Naghizadeh-Khouei & Clarke 
( 1993), the probability distribution function of the observed 
polarization angle is well approximated by a Gaussian as 
long as the Stokes parameters are Gaussian random vari- 
ables and the signal-to-noise ratio is large enough. Since the 
Faraday rotation measure between two wavelengths is pro- 
portional to the difference of the polarization angles, RM 
is also a Gaussian random function. According to our esti- 
mates, the deviations from the Gaussian distribution in RM 
do not exceed a few percent for signal-to-noise ratios exceed- 
ing unity. In Fig. ^ the signal-to-noise ratio is larger than 
3 (Beck et al. 1998). 



3 LEVEL LINES OF A RANDOM FUNCTION 

It is useful to distinguish two types of level lines of a random 
distribution of RM, and so two types of canals. Level lines 
of the first type occur at a level close to the mean value of 
RM, denoted here RM, whereas the others occur near high, 
rare peaks of RM. The two types have distinct statistical 
properties. 

Consider a random function of position F(x), with <jf 
its standard deviation, a% = {[F(x) — F] 2 }, where angular 
brackets denote averaging, F — (F(x)) the mean value of 
F(x), and C(r) its normalized autocorrelation function, 

a F C(r) — (F(x)F(x + r)} — ~F 2 , r = \r\ , 

where F(x) has been assumed to be statistically homoge- 
neous, so the autocorrelation function depends on r alone. 



3.1 Passages through a level Fo 

Now consider positions where F(x) passes through a certain 
level Fo that does not differ much from the mean value of F, 
i.e., \F — Fo\ ~ of (as we discuss below, this is the case in 
Fig. . The corresponding canals will form cells where they 
surround local extrema of Fix) or will have low curvature, 
like a contour of constant height along a valley or a mountain 
ridge. 

A useful summary of the theory of passages of a random 
function F(x) through a level Fo can be found in Sect. 33 
of Sveshnikov (19651. The theory applies to both Gaussian 
and non-Gaussian random fields, but here we only present 
results for the Gaussian case. The total probability of both 
upward and downward passages is given by 
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P(Fo) = exp 

7T OF 



(F - Ff 



2a F 



(3) 



i/2A 2 - |RM| 



where ovf is the standard deviation of the spatial deriva- 
tive of F(x). The mean number of the level lines per unit 
length is simply given by no = P(Fq). Since o%j F /o F = 
— d 2 C(r)/dr 2 | _ Q = 2/l T , where It is known as the Taylor 
microscale (e.g., Sect. 6.4 in Tennekes & Lumley 1972), this 
yields the following expression for the mean separation of 
level lines: 



xo — n — —= It exp 
v2 



(F - Ff 



2a% 



(4) 



In application to the level lines of RM discussed below, 
one should allow for passages through both the levels +Fq 
and —Fq because the canals occur wherever |RM| = RMo. 
Then Eq. ^ contains an additional term on the right-hand 
side, proportional to exp[— (Fo + F) 2 /2a F ]. This term can 
be significant if F ~ 0. 



3.2 Level lines around high peaks 

In this section we consider level lines that occur near high 
extrema of F(x) and so tend to be closed lines enclosing a 
relatively small region, like a contour of constant height near 
a high mountain peak. In the astrophysical context, statis- 
tics of high peaks have been explored in studies of structure 
formation by gravitational instability in a random density 
field (Peebles 1984; Bardeen et al. 1986). 

Consider the neighbourhood of a position where F(x) 
has a high peak, and assume that F(x) > 0. For conve- 
nience, let this position be x = 0, so that VF(0) = and 
F(0) = vgf, where v (2> 1) is the dimensionless height of 
the peak. Then F(x) at a given x near the peak is a ran- 
dom variable. The main result useful for our purposes is 
that, for a Gaussian random function, the distribution of 
F(x) is Gaussian with the mean value 



{\F(r) - F\) = F(0)C(r) 



(5) 



where r is the distance from the peak. In other words, the 
mean profile of a Gaussian random function near a high peak 
is given by its autocorrelation function. Then the distance 
of the level line from the peak, ro, follows from Eq. (J^J with 
(\F(r) - F\) = \F - F\ and C{r) = C(r ). 

Assuming that the autocorrelation function has a 
power- law form in the relevant range of scales, C(r) = 
1 — Ar^ with certain constants A and /j, (/i > 0), we ob- 
tain the following expression for the mean separation of the 
level lines, xq = 2ro: 



xo = 2 



3" 



\F - F\ 



(6) 



4 THE MEAN ANGULAR SEPARATION OF 
CANALS 

The application of the above results is straightforward. Now 
F = |RM(x)|, where x is position in the plane of the sky. 
Equations @ and © with F = |RM| and RM = nn/(2X 2 ) 
show that the angular separation of the canals depends on 
the mean value and standard deviation of RM, and on A, 
via 



\F -F\ = 

C"F CRM 

According to Eq. 0), xo decreases as A grows for A 2 < 
n7r/(2|RM|) and increases at larger A. Interestingly, Eq. © 
shows the opposite behaviour. 

To discuss the properties of the canals, we use the map 
of Fig. where the canals are very prominent and where 
we can obtain reliable estimates of RM and ctrm ■ The mean 
value of RM in the direction of M31, obtained from the ro- 
tation measures of extragalactic radio sources in that part 
of the sky, is RM fg ~ -90radm" 2 (Beck 1982; Han, Beck & 
Berkhuijsen 1998). This value of RM, produced in the Galac- 
tic foreground, can be identified with the intrinsic Faraday 
rotation measure, or the Faraday depth in that direction, 
0.81 J n c B\\ ds with n c the electron density, B\\ the line- 
of-sight magnetic field component, and L the pathlength 
through the magneto-ionic medium. However, for polarized 
emission arising within the Milky Way the maximum ob- 
servable RM would be iRMf g if polarized emission from 
all distances through the Milky Way were visible (because 
of differential Faraday rotation). The observable RM can 
be further reduced by depolarization. Then RM that enters 
Eq. Q is just a part of iRMf g . An estimate of the depth 
through the Milky Way over which the observed polarized 
emission is produced and the corresponding effective value 
of RM can be obtained as follows. 

The dominant mechanism of depolarization, apart from 
differential Faraday rotation, is internal Faraday dispersion 
due to turbulence in the ISM. The standard deviation of RM 
in the field of M31 calculated between A6cm and All cm is 
urm — lOradm -2 at a scale of 5' (Berkhuijsen et al. 2003). 
Using Eq. (34) of Sokoloff et al. H1998L we conclude that this 
mechanism reduces the degree of polarization at A20.5 cm by 
about 16%. Therefore, the effective geometric depth visible 
at this wavelength is about 0.84 of the path length through 
the magneto-ionic layer of the Milky Way in that direction. 
(Thus, the region visible in polarized emission at A20.5 cm 
in the direction of M31 is about 2kpc deep, assuming that 
the scale height of the magneto-ionic layer is 0.8 kpc.) Then 
the observable Faraday rotation measure produced within 
the Milky Way at A20.5 cm can be estimated as 

RM ~ | 0.84 RMf g ~ -38radm" 2 . 

The smallest value of RM at which complete depolariza- 
tion occurs at A = 20.5 cm is RMo ~ 37radm~ 2 from 
Eq. ©. Since RMo — |RM| is even smaller than <trm, re- 
sults of Sect. 13.1 1 apply and Eq. @ yields xo ~ 2.2Zt for the 
mean separation of canals. Our estimate of the mean sepa- 
ration of canals in Fig. Q gave x ~ 5', resulting in It ~ 2'. 
This result is discussed in Sect.|K| 

A useful diagnostic of the effects discussed here is the 
displacement of the canals as A changes so that RMo of 
Eq. (0 changes by ARM < RM . Equation Q yields the 
following estimate for the relative displacement of a canal 
(equal to half the increment in xq)- 



Ax 

x 



ARM,, ■ 



RMo - IRM! 



2cr 2 
RM 



where we have neglected terms quadratic in ARMo. For 
the observations of Gaensler et al. (2001) where |RM| ~ 
— 13radm -2 , <jrm — 45radm -2 (B. Gaensler, private com- 
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munication) and ARMo ~ 4.8 rad m~ 2 , we obtain Ax/xo ~ 

0. 024. We estimate the typical separation of canals in the 
maps of Gaensler et al. as xo — 7', which yields Aa; ~ 0.'2, 

1. e., about 1/6 of the beam width in these observations. 

Alternatively, the displacement of an individual canal 
between wavelengths Ai and A2 can be estimated if the 
gradient of RM, VRM, is known as Ax = ^7r(Aj~ 2 — 
A^~ 2 ) IVRMp 1 . For the observations of Haverkorn et al. 
(2000) , where RM varies by about 5 rad m -2 across the beam 
width of 4', we obtain Aa; ~ 0.'3, i.e., less than 10% of the 
beam width. 

For the sake of illustration, let us estimate the sepa- 
ration of canals that occur near high peaks of |RM|. Since 
RM - |RM| < o-rm in the field of M31 at A20.5 cm, the the- 
ory described in Sect. 13.21 is not useful here and we scale 
the result to All cm. Minter & Spangler (.19961 suggest the 
following form of the structure function of RM fluctuations 
in the directions around (l,b) = (144°,— 21°), i.e., not far 
from the field shown in Fig. □ D(r) ~ 5910 r 5/3 rad 2 m~ 4 
at r <0°14 (~ 8'), with r the angular separation in degrees. 
With C(r) = 1 - D(r)/(2cr RM ), Eq. © yields for Fig.0 
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^ ~ 2 



47rad 2 m~ 4 

/0.11m\ 2 
V A J 



130 rad m 

W RM 

|RM| 



130 rad m~ 



3/5 



(7) 



5 DISCUSSION 

As shown above, the mean separation of canals in Fig. 
is related to the Taylor microscale, a measure of the cur- 
vature of the autocorrelation function of RM at small sep- 
arations. In the present context, the Taylor microscale can 
be tentatively identified, e.g., with the thickness of mag- 
netic ropes, or with the size of blobs or the thickness of 
filaments in the thermal electron distribution. We empha- 
size that the canals are not directly related to any ropes or 
filaments in the magneto-ionic medium. However, the sta- 
tistical properties of the canals are sensitive to small-scale 
structures in the medium. Ropy structures of magnetic fields 
arise naturally from fluctuation dynamo action (Zeldovich, 
Ruzmaikin & Sokoloff 1990; Subramanian 1999), and fil- 
amentary structures in the gas components are abundant 
in the ISM (e.g., Koo, Heiles & Reach 1992). The result- 
ing estimate of the the rope thickness, 2', corresponds to 
about 0.6 pc at a distance of 1 kpc. Faraday rotation from 
such structures would be difficult to observe directly as they 
would need to have n e B\\ > 20 fiG cm -3 in order to produce 
|RM| ><7rm ~ 10 rad m" 2 . 

When combined with Faraday dispersion, depolariza- 
tion by differential Faraday rotation leads to a smooth de- 
pendence of p on RM, where p does not entirely vanish 
at the local minima (Sokoloff et al. 1998). This leads to 
canals where polarized intensity does not vanish but is only 
strongly reduced; still, they preserve, in an approximate 
manner, the properties discussed above. 

Canals have been observed in many directions in the 
Milky Way. In particular, Duncan et al. ( 1997J have detected 
them at the wavelength A12.6cm in the southern Galactic 
plane. As follows from Eq. (|5J, the value of RMo at that 



wavelength is at least 100 rad m -2 , and so the canals are ob- 
served because the mean Faraday rotation measure is suit- 
ably high in that direction. 

A signature of the mechanism discussed here is the dis- 
placement of the canals when the wavelength varies. Canals 
produced by discontinuities in magnetic field or Faraday ro- 
tation measure would not shift with wavelength, whereas 
those produced by differential Faraday rotation must shift 
as described in Sect. 3] However, canal displacements con- 
sistent with our model may be too small to be detected in 
the data available. 

Another interesting type of ghosts in polarized emis- 
sion is discussed by Gaensler et al. (2001). These ghosts 
arise in interferometer observations because an extended, 
polarized (and unpolarized) background, uniform in inten- 
sity (to which an interferometer is insensitive), is missing. 
However, small-scale structures in polarization angle (i.e., in 
the Stokes parameters Q and U) are not missed and cause 
variations in the observed polarized intensity but not in the 
total synchrotron intensity. This can result in blobs and fil- 
aments of enhanced or reduced polarized intensity. These 
effects may lead to extreme fractional polarization, possibly 
even above 0.75. Wieringa et al. (1993) describe polarized 
filaments invisible in total intensity, observed with the West- 
erbork interferometer. We note that ghosts of this type also 
occur in single-dish observations (Uyaniker et al. 1999a). 
Except for some fields of Uyaniker et al. 11999a! , such mea- 
surements are rarely absolutely calibrated in both total and 
polarized intensities. Therefore, a baselevel is subtracted to 
obtain a zero level for the emission. Such ghosts may be 
called 'base-level ghosts '. Furthermore, missing spacings and 
baselevel corrections have the same effect as a missing (i.e., 
depolarized) quasi-uniform, magneto-ionic layer. They lead 
to a shift in the positions of the canals and possibly to their 
(dis) appearance, without changing their statistical nature. 
As base-level ghosts can be produced at one's will in any 
data by adjusting the zero level, this opens new possibilities 
to study the magneto-ionic medium using the ghosts. 

In summary, both Faraday ghosts (canals) and base- 
level ghosts (blobs, depressions, filaments) are features in 
the observed polarized emission without a counterpart in 
the total emission. As they do not arise from any real struc- 
tures in the ISM, we call them 'ghosts'. The canals are due to 
depolarization by differential Faraday rotation in the turbu- 
lent magneto-ionic medium. They are signs of the turbulent 
structure of the ISM over a range of scales, and so represent 
a useful probe of interstellar turbulence. 
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